Electron energy spectrum and the Berry phase in graphite bilayer 
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We emphasize that there exist four Dirac-type points in the electron-energy spectrum of a graphite 
bilayer near the point K of its Brillouin zone. One of the Dirac points is at the point K, and three 
Dirac points lie nearby. Each of these three points generates the Berry phase ir, while the Dirac 
point at K gives the phase — ir. It is these four points that determine the Berry phase in the bilayer. 
If an electron orbit surrounds all these points, the Berry phase is equal to 2n. 

PACS numbers: 73.63.Bd, 71.70.Di, 73.43.Cd, 81.05.Uw 



In the recent Letter— McCann and Fal'ko considered 
the electron-energy spectrum of a graphite bilayer in the 
vicinity of the point K of its Brillouin zone and stated 
that in this bilayer the low-energy electronic excitations 
correspond to chiral quasiparticles with a parabolic dis- 
persion exhibiting Berry phase 2ir. This value of the 
Berry phase explains the unconventional Hall effect ob- 
served in the bilayer 2 . Here we refine some details of the 
spectrum and show that in reality this Berry phase in the 
bilayer is generated by four Dirac points of its spectrum. 

First of all we point out that the electron spectrum of 
the bilayer— can be easily obtained from the well-known 
Slonzewski-Weiss-McClure spectrum^ of bulk graphite if 
one puts cos£ = 0.5 and 72 = 75 = (the parame- 
ters 72, 75 describe the interaction of the atoms in the 
next-nearest-neighbor layers of graphite that are absent 
in the bilayer; £ is the dimensionless wave vector perpen- 
dicular to the graphite layers). The Slonzewski-Weiss- 
McClure model^ describes the wave-vector dependence 
of four electron energy bands of graphite E(k) in the 
vicinity of the vertical edge HKH of its Brillouin zone. 
These bands can be found from the forth-order secular 
equation: 



det H - E = 0, 



where the Hamiltonian matrix H has the form 
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Here the following notations have been used: 
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FIG. 1: Shown is the dispersion law E(k x , k y ) of the electrons 
(e) and holes (h) in the graphite bilayer near the point K. The 
spectrum contains four Dirac points. At the energy E c ~ 2 
meV all the electron cones merge. Note that without neglect- 
ing the parameter 74 ~ 0.044 eV, the spectrum E(k x ,k y ) is 
slightly asymmetric relative to the plane E = 0: The Dirac 
points in the central and side cones are at energies E — 
and E — 274717I/70 ~ 0.15-B C , respectively. Thus, in con- 
trast with graphene, the graphite bilayer is a semimetal, and 
changing the Fermi level in it, one cannot obtain the concen- 
tration of the charge carriers less than ~ 10 s cm" 2 . Scales of 
the plot are distorted for clarity. 



#33 = 7 3 re ia c, 

where a is the angle between the direction of the vector 
k and the TK direction in the Brillouin zone; T = 2 cos£; 
£ and C are dimensionless wave vectors in the direction of 
the z-axis (i.e. HKH axis) and in the basal plane, respec- 
tively: £ = (n/ 2)(k z /\KH\), ( = (2n/V3)(k ± /\TK\); 

fcj 



k^; k is measured from the point K. The 



2 



parameter 70 which describes the interaction between 
neighbor atoms in a graphite layer is sufficiently large as 
compared to the other parameters 74, A which describe 
relatively weak interactions between atoms in different 
graphite layers* 4 - As it was mentioned above, the spec- 
trum of the bilayer is obtained if one puts T = 1 and 
72 = 75 = 0. Note that in this way one can allow for the 
small parameter 74 that was neglected in Ref. [l|. 

In the interval E c <C \E\ <C 71 equations CO)-© lead 
to the approximate formula for the two low-energy bands 
of electrons (e) and holes (h) in the bilayer 



7o ,1 



7i 



(4) 



that exhibits a quadratic dependence on k± discussed in 
Ref.fl Here E c = (71/4X73/70) 2 « 2 meV, and the signs 
plus and minus correspond to the electrons and holes, 
respectively. 

At energies \E\ < E c the role of the so-called trigo- 
nal warping- increases, and this warping breaks the line 
E = constant in the k x -k y plane into one central and 
three side pockets^ see Fig. 2 in Ref. U. We emphasize 
here that each of these pockets contains a point at which 
the electron and the hole bands contact, and near all 
thesepoints the spectrum is linear in k, Fig. 1, see also 
Refs.[6.7.8. Thus, similarly to bulk graphite^ and in con- 
trast with graphene, in the graphite bilayer near the point 
K there are four points of the Dirac type. The central 
contact point coincides with K, while the three side con- 
tact points are at a distance of (V3/27r) (7371/70 )|rif| ~ 
0.005|riT| from K where \TK\ is the distance between 
the point K and the center T of the Brillouin zone. 

Before analyzing the Berry phase in the bilayer, we 
point out a general property of this phase for closed semi- 
classical electron orbits in crystals with inversion symme- 
try in the magnetic field H. The Berry phase <&b for such 
an orbit T lying in the Brillouin zone and belonging to a 
band is given b y 10 i n 
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where the direction of the integration is determined by 
the vector [H x dE°/d\t], SI is the intraband matrix ele- 
ment of the periodic (in k) part of the coordinate opera- 
tor in the crystal momentum representation^ 



n(k) 



dru£ (r)V k Uko(r), 



(6) 



and Wko( r ) is the periodic factor in the electron Bloch 
wave function of the band 0, V'ko( r ) = exp(zkr)uko(r). 
This Berry phase manifests itself in the Onzager-Lifshitz- 
Kosevich quantization condition 13 ' 14 for energy levels e 
of an electron in the magnetic field, 



S(e,k H ) 



2iTeH 

he 



(rc + 7), 



(7) 



where S is the cross-sectional area of the closed orbit in 
the k space; kn is the component of k along the magnetic 



field H; n is a large integer (n > 0); e is the absolute value 
of the electron charge, and the constant 7 is given by the 
formulas 1 ^ 

1 <I>R 

^=2-2? ' < 8 » 

If one reverses the direction of the magnetic field, the 
direction of the integration in Eq. ([5]) is also reversed, 
and the Berry phase changes its sign, while 7 changes by 
$3/77. But the electron spectrum in the magnetic field 
has to be invariant under this transformation. Since in 
the semiclassical approximation the quantity 7 is defined 
up to an integer, we conclude that the Berry phases for 
such the orbits are always multiple of irr^ This property 
of the Berry phase agrees with the results of Ref. [Ill- 
It was shown in that paper that if the electron orbit T 
surrounds band-contact lines (Dirac points in the two 
dimensional case), each of the lines (the Dirac points) 
contributes ±7r to $b ■ If the orbit does not surround the 
band-contact line (the Dirac points), the Berry phase is 
equal to zero. This property of the Berry phase also 
means that $ b can change only abruptly when the crys- 
tal potential is perturbed, and a small variation of this 
potential can, in principle, lead to an essential change of 
the Berry phase. 

McCann and Fal'ko derived the Berry phase for the 
electron orbits in the bilayer from an effective Hamil- 
tonian that leads to the parabolic spectrum How- 
ever, it is clear that parabolic spectra are idealization, 
and they never occur in crystals, and at least a small 
warping of these spectra always exists. As it was men- 
tioned above, this small warping may essentially change 
the Berry phase, and in principle, different symmetries 
of the small warping may lead to different values of the 
Berry phase since the number of the Dirac points depends 
on symmetry of a crystal. In the case of the bilayer the 
trigonal warping generates the three additional side Dirac 
points and changes the type of the central band-contact 
point. Thus, for the derivation of the Berry-phase value 
to be justified, it is necessary to consider the real sym- 
metry of the spectrum in the bilayer and to take into 
account all the band-contact points of this spectrum. 

Although each of the four Dirac point in the bilayer 
generates the Berry phase ±7r, it is necessary to find the 
signs of these phases. This can be done using the ap- 
proach of Ref. [Tg in which the effect of a small spin-orbit 
interaction on the Berry phase was investigated. It is im- 
portant that if this interaction is weak, it does not change 
the Berry phase. But the interaction enables one to fix 
the sign of the phased This is clear from the following 
considerations: Without the spin-orbit interaction the 
quantity [Vk x SI] is singular at the Dirac points. The 
interaction lifts the band degeneracy, and [Vk x SI] be- 
comes a smooth function which can be calculated using 
the formula 1 ^ [see also expressions (9), (10) in Ref. [171 
and formula (6.6) in Ref. [l8| . 
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FIG. 2: The field fl^, fc,,) (arrows) for the electron band e 
in the bilayer. The length of the arrows is proportional to |fl|. 
Small circles mark the positions of the Dirac points in the k x - 
k y plane, and the solid line shows an electron orbit with E > 
E c . Note that near the central and the side Dirac points the 
vector fi "circulates" in opposite directions. The components 
ki are given in the dimensionless units f< = (2n/\3)(ki/\TK\) 
where \TK\ is the distance between the point K and the center 
of the Brillouin zone, F. 

that is completely independent of a gauge of the elec- 
tron Bloch wave functions. Here i>q ; are iterband matrix 
elements of the velocity operator for the bands and I 
at a point k. Taking into account that Eq. (J5]) may be 
rewritten in the form 

$ B = / dk x dk y [v k x n]„ (io) 

Js 

where the integration is over the surface S spanning the 
orbit r, one can unambiguously calculate the Berry phase 
for each of the Dirac points in the bilayer. 

The effect of the spin-orbit interaction on the spectrum 
of graphene was considered in Refs. [iJSil and the fol- 
lowing Hamiltonian for the electron and hole bands in 
the vicinity of the point K was obtained 

H = v(k x a x + kytjy) + A so a z s z , (11) 

where <Ji are the Pauli matrices describing the two bands, 
v is a matrix element proportional to 70, 2A so ~ 0.2 meV 
is the gap due to the spin-orbit interaction, and s z is the 
Pauli matrix representing the electron's spin. Since the 
interaction between carbon atoms in a graphite layer is 
larger than the interaction between the atoms in different 
graphite layers (i.e. 70 S> 7i), we use Hamiltonian (|11|) 
of a single layer to include the spin-orbit interaction in 
Hamiltonian ([2]), (J3j) of the bilayer. Then, in the leading 
order in the small parameter A so /7i the matrix elements 
E3 in the third and forth lines of formula @ should be 
replaced by E 3 — A so and E 3 + A so , respectively. 

Although formulas ([9]), (fTQ|) clearly demonstrate the 
invariance of the Berry phase relative to the gauge trans- 
formations of the Bloch wave functions, it is more con- 
venient to find the field £l(k x , fe„) in the bilayer, and to 



calculate the Berry phase directly from formula (O . This 
field for a band can be obtained with the formula ; 15 > 21 

where S'(k) is the matrix reducing Hamiltonian (JSJ, © 
to the diagonal form, S + is the Hcrmitian conjugate ma- 
trix, and the subscript means that one has to consider 
the diagonal matrix element corresponding to the band 
0. Assuming the spin-orbit interaction is infinitesimal, a 
direct calculation of S leads to the field ft(k x , k y ) shown 
in Fig. 2. Note that near the central and the side Dirac 
points the vector $1 "circulates" in opposite directions, 
which means the opposite signs of the Berry phases gen- 
erated by these points. The calculation of the integral 
([5]) gives the Berry phase 7r for electron orbits surround- 
ing each of the three side Dirac points and — it for an 
orbit around the central point. Since the Berry phase of 
an electron does not depend on a size or a shape of its 
orbit in the Brillouin zone but is determined only by the 
Dirac points enclosed by the orbit ^ one finds the Berry 
phase 3 • it — tt — 2ir for electrons with energies E > E c 
(i.e., at the electron concentration N > N c 1=3 5 • 10 10 
cm~ 2 ) when their orbits surround all these four points. 
Thus, at E > E c we arrive at the same Berry phase 2n as 
in the case of the parabolic spectrum ^ in spite of the 
change in the energy-band degeneracy. The difference in 
the Berry phases for the parabolic and real spectra will 
manifest itself only in the interval —E c <E<E c &t 
low magnetic fields when many Landau levels lie in this 
interval, and hence the trigonal warping is not a small 
perturbation as compared to the magnetic energy of the 
electron. 

The above considerations also clarify existence of 
the resistivity maximum discovered in the bilayer. 2 In 
graphene the universal resistivity maximum h/Ae 2 was 
observed at zero magnetic field and low charge-carrier 
concentration N, and this maximum was explained by 
absence of localization for electrons with the Dirac-type 
spectrum^ A similar resistivity maximum was also ob- 
served in the bilayer^ and Novoselov et al£ emphasized 
that this observation is unexpected due to the parabolic 
spectrum in the bilayer. Although the parabolic spec- 
trum of Ref. [l] and the spectrum with the four Dirac 
points lead to the same Berry phase at TV > N c , existence 
of the Dirac points seems to shed a light on appearance 
of the resistivity maximum in the bilayer. The effect of 
the trigonal warping on the resistivity maximum in the 
bilayer was quantitatively analyzed in Refs. [rj|23l. 

Finally, let us point out that our result for the Berry 
phase is reminiscent of the result of Ref. d. Authors of 
that paper studied stability of Fermi points in multilayer 
graphene relative to perturbations of the crystal poten- 
tial, and they found that in the bilayer some topological 
charge Q is equal to 2 for the case of the parabolic spec- 
trum ([?]), while if one takes into account the trigonal 
warping, each of the three side Dirac points located near 
the point K has the charge 1 and the charge of the cen- 



4 



tral Dirac point at K is —1. Nevertheless, it is necessary 
to emphasize that in general case formula (11) defining 



Q in Ref. [a does not coincide with <f>s /ir [compare this 
formula with our Eqs. and (fT2|) ]. 
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